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Abstract. We review various approaches to approximating general relativistic effects in hydrodynamic simulations of stellar 
core collapse and post-bounce evolution. Different formulations of a modified Newtonian gravitational potential are presented. 
Such an effective relativistic potential can be used in an otherwise standard Newtonian hydrodynamic code. An alternative 
approximation of general relativity is the assumption of conformal flatness for the three-metric, and its extension by adding 
second post-Newtonian order terms. Using a code which evolves the coupled system of metric and fluid equations, we apply 
the various approximation methods to numerically simulate axisymmetric models for the collapse of rotating massive stellar 
cores. We compare the collapse dynamics and gravitational wave signals (which are extracted using the quadrupole formula), 
and thereby assess the quality of the individual approximation method. It is shown that while the use of an effective relativistic 
potential already poses a significant improvement compared to a genuinely Newtonian approach, the two conformal-flatness- 
based approximation methods yield even more accurate results, which are qualitatively and quantitatively very close to those 
of a fully general relativistic code even for rotating models which almost collapse to a black hole. 
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INTRODUCTION 

Gravity is the driving force which at the end of the life of massive stars overcomes the pressure forces and causes 
the collapse of the stellar core. Furthermore, the subsequent supernova explosion results from the fact that various 
processes tap the enormous amount of gravitational binding energy released during the formation of the proto-neutron 
star Therefore, gravity plays an important role during all stages of a core collapse supernova. General relativistic 
effects are important in this scenario and cannot be neglected in quantitative models because of the strong compactness 
of the proto-neutron star It is thus essential to include a proper treatment of general relativity (GR) or an appropriate 
approximation in the numerical codes that one uses to study core collapse supernovae. 

Stellar core collapse is also among the most promising sources of detectable gravitational waves. As the first-stage 
gravitational wave interferometer detectors (like TAMA300, GEO600, LIGO, or VIRGO) are already taking data, 
the interest in performing core collapse simulations has been further motivated by the necessity of obtaining reliable 
gravitational waveform predictions for detector data analysis. 

In the last two decades, a number of numerical simulations of rotating stellar core collapse were performed, paying 
particular attention to computing gravitational wave emission (see [1] and references therein). Most simulations were 
done using simple analytic equations of state (EOS), and were either restricted to Newtonian gravity or at most 
approximations of full GR. Only recently there were several successful attempts to model rotating core collapse to a 
proto-neutron star in full GR [2, 3]. Multi-dimensional simulations involving a sophisticated treatment of microphysics 
like a tabulated EOS and neutrino transport are up to now only performed in Newtonian gravity (see, e.g., [4]). 

In order to facilitate the numerical investigation of rotating stellar core collapse taking into account effects of GR 
in ways which are both accurate and comparatively easy to implement without having to resort to a fully relativistic 
formulation, in recent years a number of approximation scenarios were suggested. In this work we summarize several 
such methods and apply them to numerical simulations for modeling the collapse of rotating massive stellar cores. 
These models, which were investigated in [3] in a full GR framework, feature comparatively extreme initial masses, 
rotation states, and collapse parameters, and thus reach very high densities or even collapse to black holes. Therefore, 
this is an ideal test scenario to assess the quality of the various approximation methods of full GR. 

Throughout the article, we use geometrized units with c — G — I and Einstein's summation convention. 



HYDRODYNAMIC EQUATIONS IN FLUX-CONSERVATIVE HYPERBOLIC FORM 



In the following we present the hydrodynamic equations in flux-conservative formulation, which govern the evolution 
of matter (approximated as an perfect fluid) in a dynamic spacetime g^v We adopt the 3 + 1 formalism to foliate 
spacetime into a set of non-intersecting spacelike hypersurfaces Zf . The line element then reads 

ds^ = gijdx^dx' = -a^dt^ + Yij{dx' + li'dt){dx' + p'dt), (1) 

where a is the lapse function describing the rate of advance of time along a timelike unit vector normal to Zj, 

j3' is the spacelike shift three-vector which specifies the motion of coordinates within a surface, and Yij is the spatial 
three-metric. In the Newtonian limit g^y reduces to the flat metric g^y with the associated flat three-metric fij. 



Newtonian hydrodynamics 



As primitive (physical) hydrodynamic variables of a perfect fluid we choose the rest-mass density p, the covariant 
three-velocity v,- as measured by an Eulerian observer at rest in and the specific internal energy s. We introduce the 
foUowing set of conserved variables: 

D = p, Si = pVi, T = p{e + ^Viv'). 

Then we obtain a first-order, flux-conservative hyperbolic system of hydrodynamic equations in Newtonian gravity. 



where 7 is the determinant of 7,7. The state vector U, flux vector F\ and source vector S are given by 



(2) 



u = [D,Sj, t] , F' = [Dv',Sjv' + 5;p, (t + py] 



Here 4> is the Newtonian gravitational potential. An EOS, which relates the fluid pressure P to some thermodynami- 
cally independent quantities, e.g., P = P{p,e), closes the system of conservation equations. 



General relativistic hydrodynamics 

The hydrodynamic evolution of a perfect fluid in GR with four-velocity u^, rest-mass current = pu^, and 
energy-momentum tensor T^^ = phu^u^ -|- ''^ is determined by a system of local conservation equations, 

V^/''=0, V^7'''^=0, (3) 

where denotes the covariant derivative. The quantity h= 1+e + P/p is the specific enthalpy, and the three- 
velocity is given by v' = u'/{au^) + j3 '/ a. Along the same Unes as in the Newtonian framework, and foUowing the 
work presented in [5], we define the set of conserved variables as 

D^pW, Si=phW\, T = phW^-P-D. 

In the above expressions W = au^ is the Lorentz factor, which satisfies the relation W = 1/ \J'\'--T^vc^j. 

Then the local conservation laws (3) can again be written as a first-order, flux-conservative hyperbolic system of 
hydrodynamic equations in GR gravity. 



dt dx 
with 

U = [D, Sj, t] , F' = [Dv', Sjv' + 5jP, tv' +Pv'], S - 



Here v' =v' — P'/a, and g and 7 are the determinant of gij and Yij, respectively, with = 0!y^. In addition, F^,, 
are the Christoffel symbols associated with g^y. 



GRAVITATIONAL FIELD EQUATIONS 



Effective relativistic potential for Newtonian simulations 

Method N: Newtonian potential for a self -gravitating fluid. The Newtonian potential whose gradient acts as an 
external force in the source of the Newtonian hydrodynamic equations (2), is determined by the Poisson equation 

A<J> = 4a:p. (5) 

To approximately include effects of GR gravity, several successful attempts of increasing complexity have been 
made [6, 7, 8] to replace the regular Newtonian potential <I> (called method N here), by an effective relativistic 
gravitational potential <J>efe for a self-gravitating fluid, which mimics the deeper gravitational well of the GR case. 

Method R: TOV potential for a self-gravitating fluid. As a first approximation, we demand that <t>efr in spherical 
symmetry reproduces the solution of hydrostatic equilibrium according to the Tolman-Oppenheimer-Volkoff (TOV) 
equation. To obtain an effective relativistic potential also for a nonspherical matter distribution, we construct a 
generalized TOV potential [6] as a radial integral over angular averaged matter quantities p, e, vv, andP: 

— {m + 4Kr'¥)^. (6) 

Here v7 is the local radial velocity of the fluid. The TOV mass fn and metric function F are given by 

m{r) = 47ij^'dr'r'^p{l+e), T= ^Jl+v;^-^. (7) 

Then the effective relativistic potential Ogjf can be calculated by substituting the "spherical contribution" 

^(r) = -An[ dr'r'^ , ^ ,, (8) 
Jo \r-r'\ 

to the multi-dimensional Newtonian gravitational potential 3> by the TOV potential 3>tov: 

'^'eff = <J> - 5 -l-^TOV- (9) 

When this reference case of the effective relativistic potential, method R, is used in an otherwise Newtonian 
simulation, the collapse of a spherical stellar core to a proto-neutron star yields after core bounce and ring-down a 
new quasi-equilibrium state which is the exact (numerical) solution of the TOV structure equations of GR. 

The effective relativistic potential of method R has already been used in simulations of supernova core collapse 
in spherical symmetry with sophisticated microphysics including Boltzmann neutrino transport [6, 9] as well as for 
rotating stellar core collapse with a simple matter treatment without magnetic fields [7] and with magnetic fields [10]. 

Method A: Modiflcations of the TOV potential. In a recent comparison [9] of supernova core collapse in spherical 
symmetry it was found that the use of the effective relativistic potential in the form (9) overrates the relativistic effects, 
because in combination with Newtonian kinematics it tends to overestimate the infall velocities and to underestimate 
the flow inertia in the preshock region. As a consequence, the compactness of the proto-neutron star is overestimated. 

In order to reduce these discrepancies, several modifications of the spherical TOV potential (6) have been tested 
recently [7]. In all cases the construction of the multi-dimensional effective relativistic potential according to 
Eq. (9) remains unchanged. In the most accurate of the variations introduced in [7], method A, an additional factor F 
is added in the integrand of the equation for the TOV mass (7). Since F < 1 this reduces the gravitational TOV mass 
used in the TOV potential (6), and thus attenuates the relativistic corrections in the effective relativistic potential <I>eff. 

As a consequence, a Newtonian simulation of core collapse using the effective relativistic potential of method A not 
only reproduces the solution of the TOV structure equations for a matter configuration in equilibrium to fair accuracy, 
but contrary to method R also closely matches the results from a relativistic simulation during the dynamic phase of 
the contraction of the core. This has been demonstrated in numerical studies of collapsing stellar cores in spherical 
symmetry with sophisticated microphysics including Boltzmann neutrino transport [7]. 

The results presented in that work also demonstrate that in axisymmetric collapse calculations for rotating stellar 
cores with a wide range of initial conditions and a simple EOS, the new effective relativistic potential of method A 
reproduces the characteristics of the relativistic collapse dynamics quantitatively much better than the Newtonian 
potential for not too rapid rotation. 



Conformal-flatness-based approximations for general relativistic simulations 



In the 3 + 1 formalism, the Einstein equations split into a coupled set of 12 evolution equations for the three-metric 
Yij as well as the extrinsic curvature Kij, and 4 constraint equations. We use the ADM gauge, where 7,7 can be 
decomposed into a conformally flat term with conformal factor ^ plus a transverse traceless term: 

7ij = ^^h + hJJ, wifli hJJf'=Q, f*V,/2y/ = 0, (10) 

Method CFC: Conformal flatness condition for the three-metric. In spherical symmetry h\J = 0, i.e. the three- 
metric is conformally flat. Thus a reasonable approximation of 7,7 for scenarios which do not deviate to strongly from 
sphericity is to impose a vanishing /iJJ in Eq. (10): 



Yij = ^%j- (11) 



This is the conformal flatness condition (CFC) or Isenberg- Wilson-Mathews approximation [11, 12]. Using CFC 
and assuming the maximal sUcing condition, for which the trace of the extrinsic curvature vanishes, the 3 -M metric 
equations reduce to a set of five coupled elliptic (Poisson-like) equations for the metric components, 

(K- K'-^ \ 
phW^-P+ '^^^ j, (12) 

A(a0) = 27ta<j)^ ph(3W^ - 2) +5P-I- ) , (13) 

Ion I 



AJ3' = l6na^''S'+2(i>'''K'}S/j |^^J - (14) 

where V and A = 7'^ V, V, are the flat space Nabla and Laplace operators. Additionally, the expression for the extrinsic 
curvature becomes time-independent and reads A'iy = 2^(V, j3; + V, j3, — |7,jV^j3*^). As the CFC metric equations (12- 
14) do not contain explicit time derivatives, the metric is calculated by a fully constrained approach. 

The accuracy of the CFC approximation has been tested in various works, e.g., for stellar core collapse and 
equilibrium models of neutron stars [2, 13, 14, 15, 16, 17], as well as for binary neutron star merger [18, 19]. The 
spacetime of rapidly (uniformly or differentially) rotating neutron star models is very well approximated by the CFC 
metric (11). Its accuracy degrades only in extreme cases, such as rapidly rotating massive neutron stars or black holes. 



Method CFC+: Inclusion of the second post-Newtonian deviation from isotropy. To improve the accuracy of CFC, 
we have developed the CFC+ approximation. It consists in adding to the CFC metric (11) the second post-Newtonian 
(2PN) deviation from isotropy, 

y.j = ^% + h^^ (15) 

which includes a traceless and transverse term hj™. In the decomposition (10) hj™ is identical to up to 2PN 
order. As described in detail in [20] this correction is the solution of the tensor Poisson equation [21] 

A/i2™ = 7™ (^-4%UViU-l6n{StSj) /D*) , (16) 

where {/ is a Newtonian-hke potential which is the solution of the Poisson equation At/ = —AnD*, D* = 7/ fD, and 
'^i* = VrlfSi- The integro-differential operator ffj^''' is the traceless and transverse projector defined in the Appendix 
of [20]. Of the CFC metric equations (12-14), only the equation (13) for the lapse function a is affected by the 
corrections hjj^: 

A(a(?))= [A(a(/))],2PN^o-f*7''/!2™v,V,t/. (17) 

However, the CFC+ corrections couple implicitly to the other metric equations as their source terms depend on a. 

Calculating hjj^ by inversion of the tensor Poisson equation (16) is simplified considerably by introducing poten- 
tials that are solutions of scalar/vector/tensor Poisson equations [20]. 

Applying method CFC+ to axisymmetric simulations of pulsations in rotating neutron stars and rotational stellar 
core collapse to a proto-neutron star, it was recently found that relevant quantities show only minute differences 
with respect to method CFC [20]. By contrast, in scenarios where the CFC approximation is expected to become 
increasingly inaccurate, CFC+ is expected to yield results which are closer to full GR. 



SIMULATIONS OF THE COLLAPSE OF ROTATING MASSIVE STELLAR CORES 



Setup of models and numerical methods 

We follow the work in [3J to set up rapidly rotating initial models of massive stellar cores with baryonic rest 
masses Mr> 1.5 Mq , supported by electron degeneracy pressure and a possible radiation pressure contribution. Using 
Hachisu's self-consistent field method [22], we construct the initial model as a rotating perfect fluid in equilibrium 
obeying a poly tropic EOS, P = Kqp'>'° and e =F/[p(7o — 1)], where 70 =4/3 is the adiabatic index and the polytropic 
constant assumes values Kq = {5.0, 7.0, 8.0} x lO''* (in cgs units) with increasing contribution from radiation 
pressure. The initial core's central density is Pc = lO'*^ g cm"^. The rotation law is set by j = A^{Q.c — H), where j is 
the specific angular momentum, is the value of the angular velocity Q. at the center, and the constant A determines 
the degree of differential rotation. We parameterize the rotational state by two quantities, the rotation rate j3o which is 
the ratio of rotational kinetic energy to gravitational binding energy, and the parameter A =A/rg = {0°, 0.25, 0.1}, 
where is the stellar equatorial radius, which decreases as the differentiaUty of the rotation profile increases (with 
A=°o corresponding to uniform rotation). To initiate the collapse we reduce the pressure, P = K()p'^^{Yi — 1)/ (70 — 1), 
with lowering the adiabatic index to 71 = 1.31 < 70 and setting the polytropic constant to Ki = 5.0 x lO''*. Note that 
this procedure is different from the collapse initiation for the models presented in [2, 14, 15, 20, 23], where both the 
pressure and the internal energy are reduced. 

As in [2, 3, 14, 15, 23], during the collapse phase we use a simple hybrid ideal gas EOS [24] that consists of 
a polytropic contribution Fp describing the degenerate electron pressure and at supranuclear densities the pressure 
due to repulsive nuclear forces, and a thermal contribution f th which accounts for the heating of the matter by shocks, 
P = Pp + Ptii, where fp =Ki2P^^'^ andfth = peth(7th — !)■ To approximate the stiffening of the EOS at densities larger 
than nuclear matter density Pnuc = 2.0 x lO'"* g cm~^, the adiabatic index increases to Ji = 2.5 > Ji ■ For the adiabatic 
index of the thermal contribution we choose 7th = 1.3. The thermal internal energy is calculated as fith = £ — Cp, 
while the polytropic specific internal energy gp is determined from Pp by the ideal gas relation in combination with 
continuity conditions, which also determine the value for the polytropic constant K2 at p > p^ic (for more details, 
see [14, 24]). Of the 27 axisymmetric collapse models of rotating massive stellar cores investigated in [3], here we 
select 8 representative models. Table 1 summarizes their parameters and establishes the nomenclature of model names. 

Our evolution code, which was introduced in [14, 20], utilizes Eulerian spherical polar coordinates {r, 0} restricted 
to axisymmetry. We also assume symmetry with respect to the equatorial plane. The finite difference grid consists of 
Hr = 250 logarithmically spaced radial and ng = 45 equidistantly spaced angular grid points with a central resolution 
Arc — 250 m. A small part of the grid covers an artificial low-density atmosphere extending beyond the stellar surface. 

The hydrodynamic solver performs the time integration of the system of conservation equations, given either by 
Eq. (2) for methods N, R, and A, or by Eq. (4) for methods CFC and CFC+. This is done by means of a high-resolution 
shock-capturing (HRSC) scheme with third-order accurate PPM reconstruction (for a review of such methods in 
numerical GR, see [25]). The numerical fluxes are computed by Marquina's approximate flux formula [26]. The time 
update of the state vector U is done using the method of lines in combination with a Runge-Kutta scheme with second 
order accuracy in time. Once U is updated in time, the primitive variables are recovered either directly (for methods N, 
R, and A) or through an iterative Newton-Raphson method (for methods CFC and CFC+). 

When using methods N, R, or A we solve the multi-dimensional linear Poisson equation (5) for the Newtonian 
potential «& by transforming it into an integral equation using a Green's function. The volume integral is numerically 
evaluated by expanding the source term into a series of radial functions and associated Legendre polynomials, which 
we cut at order Z = 10. It is straightforward to solve the integral equations (6-8) of the relativistic corrections needed 
to finally obtain the effective relativistic potential Ogff . 

The metric solver for obtaining the CFC and CFCh- spacetime metric used in the GR simulations exploits the fact 
that Eqs. (12- 14), supplemented by Eqs. (16, 17) in the case of method CFC+, are written in the form of a system of 



TABLE 1. Summary of the model parameters for rotating massive stellar cores. 



Model 


M5al 


M5c2 


M7a4 


M7bl 


M7c3 


M8al 


M8c2 


M8c4 


Ko [1014] 


5.0 


5.0 


7.0 


7.0 


7.0 


8.0 


8.0 


8.0 


A 


00 


0.10 




0.25 


0.10 




0.25 


0.10 


jSo [%] 


0.89 


1.24 


3.67 


2.18 


1.27 


0.88 


2.19 


1.27 



nonlinear coupled equations for a vector of unknowns u with a Laplace operator on the left hand side, Am = S{u). A 
common method to solve such equations is to keep S(h) fixed, start with an initial guess and then solve each of the 
resulting decoupled linear Poisson equations Am = S(m*) in an iteration cycle with steps s until convergence. To 
obtain the solution of the now uncoupled, linear Poisson equations, we again employ the expansion described above. 
More details are presented in the description of Solver 2 in [15]. 

The calculation of the CFC or CFC+ metric is computationally very expensive. Hence, for these methods the metric 
is updated only once every 100/10/50 hydrodynamic time steps before/during/after core bounce, and extrapolated in 
between, as described in [14]. We also note that we have performed grid resolution tests to ascertain that the grid 
resolution specified above is appropriate. When using the CFC/CFC+ metric equations, the degrees of freedom repre- 
senting gravitational waves are removed from the spacetime. Therefore, in all simulations gravitational waveforms are 
obtained in a post-processing step with the help of the quadrupole formula. 



Results and discussion 

Comparing typical quantities that describe the collapse dynamics (like the maximum density Pmaxb, maximum 
rotation rate jSmaxb. or minimal value of the lapse ttnunb during core bounce) from our approximate GR simulations 
with those obtained in full GR [3], we find excellent agreement for both method CFC and CFC+. In the case of models 
which do not collapse to a black hole, our results for Pmaxb (which we present in Table 2) deviate from the full GR 
ones by typically less than about 2% for models with moderate core mass and moderate or rapid rotation. The deviation 
for 

Pmaxb increases to at most 20% for the most extreme model M8c4, which has a baryonic rest mass of Mr = 3.05 
and a slow rotation rate. Consequently, this model reaches ttmin = 0.29 and thus almost collapses to a black hole. 

When the less accurate approximation methods A and R are used, the values for Pmaxb as expected deviate more 
strongly from full GR and approach the ones obtained using a genuinely Newtonian simulation. Concerning the quality 
of our approximations, we find equivalent results for the values of Pmaxb> J3maxb. and Umin b (with the latter obviously 
only defined for methods CFC and CFC+). Note also that though method A compared to method R yields results closer 
to full GR in spherical symmetry, the influence of rotation reverses this behavior, as already discussed in [7]. 

A very important property of any approximation method for GR simulations is that it should correctly reproduce the 
collapse type of a specific model for a range in parameter space as broad as possible. For the investigated models, 
following [3] we denote by BH, NS, O-A, and 0-B the collapse to a black hole, a neutron star, an oscillating 
star with Pmaxb > Pnuo and an oscillating star with Pmaxb < Pnuc as end state, respectively. As shown in Table 3, 
methods CFC and CFC+ correctly capture the collapse type for all investigated cases. The identification of the end 
state of model M7a4 (as well as of model M7a3) as a black hole in [3] is most probably due to insufficient resolution 
of the inner parts of the numerical grid, which results in fall-back to the already formed proto-neutron star leading to 
its recoUapse to a black hole. We find that behavior in a resolution study of these two models. This is also supported 
by independent simulations in full GR without symmetry assumptions [27]. We thus classify model M7a4 as NS. 

Fig. 1 displays the time evolution of the maximum density Pmax and the gravitational wave signal amplitude Ajg 
(i.e. the gravitational radiation waveform; for a definition of Afg, see, e.g., [23]) for models M5al and M7a4, which 
are aU of coUapse type NS in a full GR simulation. As the results of method CFC and CFC+ are almost identical, 
we do not show the latter in Fig. 1. We emphasize that even for these rather extreme stellar core collapse models, the 

TABLE 2. Maximum density Pnmx b in ™its of 10^^ g cm"^ during core bounce for the investigated rotating 
core collapse models from simulations using methods N, R, A, CFC, and CFC+. The relative deviations 

Pmaxb/Pmaxb ^om the value obtained in a full GR simulation (given in the last column) are shown in 
parentheses. The dash denotes models for which the approximate or full GR simulation produces a black hole. 



Model I Method N I Method R I Method A I Method CFC I Method CFC+ I GR 



M5al 


3.8 


(0.58) 


6.1 


(0.92) 


5.6 


(0.85) 


6.6 


(1.00) 


6.6 


(1.00) 


6.6 


M5c2 


1.1 


(0.22) 


3.3 


(0.66) 
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(0.58) 
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(0.98) 
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(0.98) 
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M7a4 
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14 
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M7bl 


0.10 


(0.13) 


0.40 


(0.51) 


0.31 


(0.39) 


0.83 


(1.05) 


0.85 


(1.08) 


0.79 


M7c3 


1.2 


(0.13) 


5.3 


(0.58) 


4.2 


(0.46) 


9.2 


(1.00) 


9.2 


(1.00) 


9.2 


M8al 


4.5 








17 














M8c2 


0.19 


(0.04) 


1.5 


(0.28) 


9.0 


(0.17) 


5.3 


(0.98) 


5.2 


(0.96) 


5.4 


M8c4 


1.2 


(0.08) 


7.1 


(0.47) 


5.1 


(0.34) 


17 


(1.13) 


12 


(0.80) 


15 



2PN deviation from isotropy in method CFC+ remains very small, h^J ~ 0.001 -0.01 <C 1. Only for models like 
M8c4 which reach extremely high densities and low values of the lapse function, the components of Ir^J^ increase 
to 0.02. This further supports the observation that method CFC already is an excellent approximation of full GR. For 
all models the purely Newtonian simulation, method N, yields maximum densities far below the values in CFC both 
during and after core bounce. These results for rotating massive stellar cores are in accordance with the findings for 
rotating regular stellar cores in [23, 28]. There the typically higher waveform amplitudes at core bounce in method N 
compared to method CFC for many models of type NS, observed here as well in Fig. 1, are also found and explained. 

In model M5al (top panel) the post-bounce value of p^ax for approximation method A matches best the result 
with method CFC, which (with method CFC+) is closest to full GR, while method R gives too large values. A similar 
behavior was found and discussed for the rotating core collapse model A1B3G1 in [7]. The inadequacy of method R for 



TABLE 3. Collapse type of the investigated rotating core collapse models from simulations using 
methods N, R, A, CFC, and CFC+. The simulations which yield the correct collapse type obtained in a 
full GR simulation (given in the last column) are marked by a box. See the main text for the definition 
of the collapse types. For model M7a4 the full GR simulation in [3] predicts BH instead of NS as final 
state, which is apparently due to lack of resolution. 
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FIGURE 1. Time evolution of the maximum density p^ax (left panels) and the gravitational wave signal amplitude Ajg (right 
panels) for the rotational core collapse models M5al (top panels) and M7a4 (lower panels) from simulations using method N (blue 
dash-dotted lines), method R (green dotted lines), method A (red dashed lines), and method CFC (black solid lines). The thin black 
solid line in the left bottom panel shows Pmax for model M7a4 using method CFC in low resolution. 



reproducing the collapse type of certain models is evident for model M7a4 (bottom panel), where the collapse proceeds 
unhalted and thus the prompt formation of a black hole is predicted, while the other approximation methods capture 
the collapse type right. Beyond employing a good approximation of GR, it is also crucial to use a numerical grid with 
sufficiently accuracy. On a grid with a low central and global resolution of rir = 100, ne = 20, and Arc — 500 m, the 
core of model M7a4 (and many other models as well) simulated in method CFC suffers from matter fall-back after 
core bounce, recoUapses, and predicts the delayed formation of a black hole as in [27] (see left bottom panel of Fig. 1). 

SUMMARY AND OUTLOOK 

We have presented several recently introduced methods for approximating GR effects in hydrodynamic simulations of 
stellar core collapse. When applied to axisymmetric test models for the collapse of rotating massive stellar cores which 
reach high densities and rotation rates, these approximations as expected offer much better results than a genuinely 
Newtonian treatment of gravity. While the accuracy of methods using an effective relativistic potential in an otherwise 
Newtonian code degrades at rapid rotation (which can be remedied by a new formulation [8]) and these methods 
reach their Umit for high densities and/or strong rotation, approximations based on conformal flatness prove to yield 
exceUent matching with fuU GR for aU investigated models. 
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